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Abstract. In the present work we report the band structure calculations for 
the high temperature superconductor Nd2-xCe a; Cu04 in the regime of strong 
electronic correlations within an LDA+GTB method, which combines the local density 
approximation (LDA) and the generalized tight-binding method (GTB). The two 
mechanisms of band structure doping dependence were taken into account. Namely, the 
one-electron mechanism provided by the doping dependence of the crystal structure, 
and the many-body mechanism provided by the strong rcnormalization of the fcrmionic 
quasiparticles due to the large on-site Coulomb repulsion. We have shown that 
in the antifcrromagnetic and in the strongly correlated paramagnetic phases of the 
undcrdoped cuprates the main contribution to the doping evolution of band structure 
and Fermi surface comes from the many-body mechanism. 
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1. Introduction 

It is well known that the hole doping of La2_ x Sr x Cu04 and the electron doping 
of Nd 2 _ x Ce a; Cu04 shift them into the superconducting state. With the increase 
of doping concentration x the band structure undergoes dramatic changes from an 
antiferromagnetic (AFM) insulator to a normal paramagnetic metal. It is well 
established that the strong electronic correlations play one of the main roles in the 
formation of the electronic structure of High-T c cuprates. They are important especially 
for small x values and should be taken into account in the calculations explicitly. 

In the present paper to study the band structure of Nd2_ a; Ce 3 ;Cu04 we apply the 
LDA+GTB computational scheme [8]. This scheme combines ab-initio calculations 
within the Local Density Approximation (LDA) and the many-body Generalized Tight- 
Binding (GTB) approach [9]. Earlier, within the GTB method a specific many-body 
mechanism of band structure evolution for the doped Mott-Hubbard insulator was found. 
This mechanism is caused by the changes of occupation factors of many-body electronic 
terms by the quasiparticle excitations [9]. For Nd2_ x Ce x Cu04 such terms are d 10 p 6 
term (zero holes per unit cell) and hybridized d 9 p 6 + d 10 p 5 term (one hole per unit cell). 
Here we highlight also the so-called "one-electron" mechanism of band structure doping 
dependence. It originates from the changes of a crystal structure (lattice parameters, 
atomic positions) upon doping. As a result, matrix elements of d — d, p — d, and p — p 
hybridization (hoppings, and all the rest depending on the interatomic distance) vary 
with doping. In case of absence of a strong electronic correlations this mechanism will be 
responsible for evolution of band structure with doping in the framework of the standard 
band theory. That is why we call it one-electron mechanism. 

Previously, the band structure of Nd2-a:Ce x Cu04 was considered within the LDA 
calculations [TJ [2j [3] and the tight-binding approach [I]. Note, our approach is 
significantly different from these two and their simple combination. From the LDA band 
structure we extract the doping-dependent tight-binding parameters using projection 
procedure, not fitting. Then, we use a many-body tight-binding computational scheme 
called GTB approach. One should not mix up the one-electron tight-binding approach 
and many-body GTB method. The details of latter will be given in the next Section. 

To our knowledge there is no band structure calculations for the High-T c cuprates 
that take both one-electron and many-body mechanisms into account. In the present 
paper we report the results of such calculations for Nd2_ a; Ce 3 ;Cu04 and compare them 
with the simplified calculation which do not contain the one-electron mechanism. It 
was found that in the antiferromagnetic phase for small doping concentrations, x < 0.1, 
the one-electron mechanism results in a slight shift of the bottom of the conduction 
band and, simultaneously, in a small shift of the top of the valence band, thus 
retaining the value of the charge-transfer gap. For higher doping concentrations the 
paramagnetic spin-liquid phase was considered in an approximation taking the static 
spin-spin correlation functions into account. In this phase, the one-electron mechanism 
provides small changes to the bandwidth. However, the Fermi surface and critical 
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concentrations at which the Fermi surface topology changes, remain the same as in the 
case when the one-electron mechanism is disregarded. 

2. Brief description of a hybrid LDA+GTB computational scheme 

Ab-initio electronic structure calculations within the density functional theory have their 
development within the LDA approximation. This approximation does not take strong 
electronic correlations into account properly. That is why the true band structure of 
Mott insulators can not be described within LDA. We employ the LDA to calculate 
the non-interacting part of the multi-band p — d model Hamiltonian. Then, the strong 
electronic correlations enter in the framework of the GTB approach [8]. 
LDA+GTB method consists of the following steps: 

(i) Ab-initio LDA band structure calculation, finding of Bloch functions; 

(ii) Construction of the Wannier functions for the physically relevant states; 

(iii) Construction of the multi-band p — d model Hamiltonian with the parameters 
obtained from the previous two steps; 

(iv) Splitting of the multi-band p — d model Hamiltonian into a sum of inter- and intra- 
cell components and exact diagonalization of the intra-cell part in order to construct 
the many-body molecular orbitals for the unit cell; 

(v) Construction of the intra-cell Hubbard AT-operators on the basis of these molecular 
orbitals and rewriting the full Hamiltonian for the crystal in the ^-representation; 

(vi) Calculation of a quasiparticle band structure for Hubbard fermions in the framework 
of the perturbation theory with the small inter-cell hopping and interactions. 

Note, the first version of the GTB method 0, [10] contained a large number of 
unknown model parameters, which were extracted by the comparison with experimental 
data. In the generalized LDA+GTB method, all parameters of the theory are calculated 
explicitly. For a given doping concentration x we calculate the ab-initio band structure. 
One-electron mechanism of band structure doping dependence is determined by the 
dependence of the matrix elements of the inter-atomic hopping and one-electron energies 
on doping, which is due to the change of the lattice parameters with x. Many-body 
mechanism arises from the doping dependence of the occupation factors. Thus, one- 
electron mechanism takes place for the ordinary tight-binding method while the many- 
body mechanism appears within the GTB method as the effect of the strong electronic 
correlations. 

3. Electronic structure and model parameters of Nd 2 _ x Ce x Cu04: LDA 
results 

Nd2Cu04 crystallizes in tetragonal structure with the symmetry space group I4/mmm 
also called the T'-structure. Lattice parameters are a = b = 3.94362A, and 
c = 12.1584 A[TT]. Cu ions in positions 2a, (0,0,0), are surrounded by four ions 
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Figure 1. The comparison of the electronic dispersion obtained within the LDA 
calculations (gray dashed curves) and the electronic dispersion of the effective non- 
interacting 3-band Hamiltonian for the NMTO orbitals basis set (black solid curves) 
for different Ce concentrations a;. 

of oxygen 01, which occupy positions 4c, (0,0.5,0). Ions of Nd in positions 4e, 
(0,0,0.35112), have eight nearest neighbors of oxygen 02 in positions 4d, (0,0.5,0.25). 
In comparison with the high-temperature tetragonal structure of La2Cu04, in the T'- 
structure of N0I2CUO4 the apical oxygen ions around Cu ions are absent. With Ce 
doping the symmetry group of Nd 2 „ x Ce x Cu04 remains the same whereas the lattice 
parameters and z-coordinate of Nd positions are changing (see Tabled]) [TT1 [T2] . 

Electronic structure of Nd2_ x Ce a: Cu04 at < x < 0.3 was calculated using the local 
density approximation. To this end, the linearized muffin-tin orbitals (LMTO) method 
in the tight-binding approach within atomic spheres approximation (TB-LMTO-ASA) 
was applied [131 E] • The 4/ states of Nd were considered as semi-core states, because 
they are well localized and are situated far below the Cu-d states [15]. The LDA band 
structure of Nd2_ x Cea;Cu04 along the high symmetry directions of the Brillouin zone is 
shown in Fig. [1] with gray dashed curves. Bands formed by the hybridized 3d copper 
and 2p oxygen states has a width of approximately 9 eV. As a result of hybridization 
between d x 2_ y 2 copper orbital and appropriate p x>y orbitals of plain oxygen (01) the 
bonding bands are located at energies —5... — 6 eV, while the antibonding bands crosse 
the Fermi level. These hybrid orbitals determine the non-interacting Hamiltonian of the 
so-called 3-band p — d model [16], Q2] . 

To calculate hopping integrals for the 3-band p — d model the NMTO (muffin-tin 
orbitals of N-th order) [18] method was used. For the three physically relevant hybrid 
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Table 1. Crystal 


structure 


parameters (A) for 


Nd 2 _xCe i 


; CuC>4 for 


different 


concentrations (see 


Refs. (TTJQ3). 










Lattice constant 


x = 0.00 


0.05 


0.10 


0.15 


0.20 


0.30 


a 


3.94362 


3.94056 


3.94071 


3.94224 


3.94295 


3.94288 


c 


12.1584 


12.1130 


12.0945 


12.0603 


12.030 


12.0288 


z(Nd) 


0.35112 


0.3519 


0.3523 


0.3527 


0.3531 


0.3531 



Table 2. Hopping integrals and one-electron energies (eV) as a function of Ce 
concentration for Nd2_ 2; Ce a ;Cu04 obtained by the NMTO method. Here x 2 , p x , p y 
denote the Cu-d x 2 _ y 2 and 0\-p x<y orbitals respectively. 







x = 0.00 


0.05 


0.10 


0.15 


0.20 


0.30 


Energy 


E{x 2 ) 




-2.2855 


-2.2847 


-2.1760 


-2.4215 


-2.3507 


-2.3234 


E( Px ) 




-3.2935 


-3.3064 


-3.2829 


-3.2607 


-3.2800 


-3.2957 


Hopping 


Direction 














t(x 2 ,p x ) 


(0.5, 0) 


1.1216 


1.1454 


1.1665 


1.1614 


1.1726 


1.1719 


t'(x 2 ,p x ) 


(0.5, 1) 


-0.0504 


-0.0359 


-0.0211 


-0.0202 


-0.0166 


-0.0201 


t"(x 2 , Px ) 


(1.5, 0) 


0.0834 


0.0921 


0.1173 


0.1130 


0.1203 


0.1126 


t'"{x 2 ,p x ) 


(1.5, 1) 


-0.0149 


-0.0083 


0.0015 


0.0090 


0.0153 


0.0099 


t(Px,Py) 


(0.5, 0.5) 


0.8320 


0.8389 


0.8381 


0.8365 


0.8386 


0.8304 


t'(Px,Py) 


(1.5, 0.5) 


0.0266 


0.0331 


0.0452 


0.0450 


0.0469 


0.0388 



bands (Cu-d x 2_ y 2 and 01-p X}V ) the NMTO orbitals were constructed. Corresponding 
band dispersions are shown in Fig. [1] by solid black curves. Observe, that NMTO 
orbitals are almost perfectly reproduce the LDA calculated bonding and antibonding 
bands formed by a hybrid d x 2_ y 2 copper and p XtV oxygen orbitals. Thus, within this 
NMTO basis we have obtained an effective few-orbital Hamiltonian, the 3-band p — d 
model, which we were looking for. 

Note, we do not take into account the Cu-4s orbitals explicitly. Their importance 
for the bilayer cuprates was shown in Refs. [Hl[6j[7]. For bilayer structures, the 4s states 
connect two neighboring Cu0 2 -layers. It is not the case for Nd2_ x Ce x Cu0 4 , the system 
where there is only one Cu02-layer per unit cell. Moreover, due to the absence of apical 
oxygens this layer becomes even more two-dimensional. There is also a contribution 
from the Cu-4s states to the in-plane hopping integrals pE]. However, in paper [5] it was 
shown that one can do the so-called Loewdin transformation and end up with the 3-band 
p — d model with the renormalized parameters. Since our tight-binding parameters are 
not fitted but calculated within the NMTO method from the full-band LDA results, 
such renormalization is automatically performed. Thus, we partially take into account 
the Cu-4s states implicitly. 

Having applied the Fourier transform to this effective Hamiltonian in momentum 
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space we obtain the real space hopping integrals depending on the distances between 
atoms (see Table [2]). From Fig. [1] and Table [2] one can conclude that hopping integrals 
change with doping slightly. It allows us to assume that the one-electron contribution to 
the evolution of Nd2_ x Ce x Cu04 electronic structure with increase of Ce concentration 
is not substantial. 

For the 3-band p — d model we also need the values of Coulomb repulsion U and 
Hund's exchange parameter Jh for Cu ions. In the paper [19] they were calculated for 
copper in the L^CuC^ compound by a super-cell method [20] and equal to: U = 10 
eV, Jh = 1 eV. We presume these values to be doping independent and use them in the 
present paper for the Nd2_ x Ce x Cu04 compound. 



4. LDA+GTB results in the AFM phase 

For the wide range of doping concentrations Nd2- x Ce x Cu04 remains in the AFM phase. 
Therefore, we will consider the evolution of the band structure with doping in the AFM 
phase first. Within the GTB method for this phase we use the Hubbard-I approximation 
[2Tj . though the diagram technique for the A-operators [22], [23], [21] allows one to go 
beyond this approximation. 

To write down the model Hamiltonian we use the Hubbard A-operators |25j : 
Ay <->• X™' n = \n) (n'\. Here index a <-> (n,n') enumerates quasiparticles with 
energy uj a = e n (N + 1) — e n >(N), where e n is the n-th energy level of the A-electron 
system. The commutation relations between A-operators are quite complicated, i.e. 
two operators commute on another operator, not a c-number. Nevertheless, depending 
on the difference of the number of fermions in states n and n' it is possible to define 
a quasi-Fermi and a quasi-Bose type operators in terms of obeyed statistics. There 
is a simple correspondence between A-operators and the single-electron annihilation 
operators: ar\ a = J2lXtr(oi)Xf, where the coefficients 7A CT (a) determine the partial 

a 

weight of the quasiparticle a with spin a and orbital index A. These coefficients are 
calculated straightforwardly within the GTB scheme. 

In the Hubbard-I approximation the dispersion equation for the band structure of 
the Hubbard fermions in the AFM phase with the sublattices A and B is the following 
HQ]: 



o, (i) 

AA' 

where is the intra-cell local energy of the Hubbard's fermion, T^y(k) is the Fourier 
transform of the matrix element of the intra-cell hopping between the one-electron 
orbitals A and A'. The occupation factor, F^, is equal to the sum of the occupancies 
of initial and final many-body states, the transition between which is described by the 
operator XJ. 

Since the structure of dispersion equation (pQ) is similar to the equation in the 
one-electron tight-binding approach, our method was called Generalized Tight-Binding 




Figure 2. Quasiparticle's energy momentum dependence calculated within the 
LDA+GTB method for x = (a), x = 0.05 (b), x = 0.1 (c), and x = 0.15 (d). 
The dotted curve corresponds to the calculation without one-electron mechanism of 
concentration dependence. The solid curve represents results of calculations in which 
both many-body and one-electron mechanisms were considered. 



(GTB) method. At the same time, there are important distinctions of equation (pQ) from 
its one-electron analogue. Namely, the local energies are calculated with explicit 
consideration of strong electronic correlations inside the cell. Thus, corresponding 
occupation factor can have non-integer values and it depends on temperature and doping 
concentration. As a result, the doping dependence of the electronic structure is not 
described by the rigid band model, so the effect of doping is not only due to the shift 
of a chemical potential for a given band structure. There are bands with the spectral 
weight proportional to x for small doping concentrations. These bands appear inside 
the dielectric gap near of the bottom of the conduction band for n-type cuprates and 
near of the top of the valence band in hole-doped cuprates. These states are called "in- 
gap states" 0, [101 [26] . Note, they are not an impurity states of a doped semiconductor 
formed in the presence of a defect since no defects are present in our model. An example 
of the in-gap states near of the bottom of the conduction band in Nd 2 _ x Ce^CuC^ will 
be discussed later in this section. 

The results of the LDA+GTB calculations of the dispersion using equation (JTJ) in 
the Neel state are shown in Fig. [2j In each figure for a given doping concentration the 
results of two calculations are shown: (i) with the model parameters calculated for x = 
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(i.e. in absence of one-electron mechanism but with the many-body mechanism taken 
into account), and (ii) for the model parameters which depend on the concentration 
(see Table E]), also taking the many-body mechanism into account. In the former case, 
the doping dependence is determined only by the doping-dependence of the occupation 
factors, F^, while in the latter case the dependence of matrix elements T^y(k) on 
doping was also taken into account. 

The results of the electronic structure calculations for an undoped compound, 
Fig.[2]^a), reproduce the main effects of the strong electronic correlations in this material. 
On the bottom of the conduction band and on the top of the valence band there are 
in-gap states with the spectral weight proportional to the doping level [HI QUI EE] . Upon 
increase of doping concentration the in-gap state at the bottom of the conduction band 
becomes dispersive with non-zero spectral weight [see Fig. [2](b-d)]. 

For each concentration one can notice that "switching-off" the one-electron 
mechanism leads to a shift of the top of the valence band and the bottom of the 
conduction band. This shift is almost uniform and its value is very small. Also, this 
does not prevent appearance of the in-gap states. Thus, one can conclude that such fine 
tuning of the Hamiltonian parameters gives just a shift of the electronic structure as a 
whole in the vicinity of the dielectric gap. This is most probably not very important for 
the physics of the High-T c materials. 



5. Paramagnetic phase 

To treat a spin-liquid phase, the multiband p—d model Hamiltonian was mapped onto an 
effective low-energy model [8] . Parameters of this effective model are obtained directly 
from the ab-initio parameters of the multiband model. The low-energy model for n-type 
cuprates is the t — t' — t" — J* model (t — f — t" — J model with the three-site correlated 
hoppings) with the following Hamiltonian in the hole representation: 

H = J2(e -»)X^+ £ t fg XfX /+J2J fg (s f S g -\n f n 9 )+Hz.(2) 

Here \i is the chemical potential, Sf is the spin operator, 5^" = Xj ,a , Sj = X"J ,a , 
Sj — ~ (Xj' a — Xj' a ^j, nf = X a f u is the number of particles operator, Jf g = 2t'j g /E ct 

is the exchange parameter, E ct = 2 eV is the charge-transfer gap. In the notations 
of Ref. [8] the hopping matrix elements tf g corresponds to — t°j° g , and tf g = —t°f g . 
Hamiltonian H% contains the three-site interaction terms: 

n i — 2^ V A / A m A g ~ A f A m A g ) ■ K 6 ) 

In the considered case there is only one Fermi-type quasiparticle, a = (0, a), with 
l\a{ct) = 1, and the Hamiltonian in the general form in the momentum representation 
is given by: 

ff = E(^-^)^r+EEft^+E E v p f™'xfx;4x§.(4) 

k,a k a >P ".ft ", "' 
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The Fourier transform of the two-time retarded Green function in the energy 



representation, G\(k,E) 



k\a 



Green function, 



D(k,E) 



a/3 



k 



k\a I I E 

/3t 



xt 



can be rewritten in terms of the matrix 



(5) 



G x (k,E) = Y.lxAuhlMD a P(k,E). 

a,(3 

The generalized Dyson equation for the Hubbard X-operators [21] in the 
paramagnetic phase {{Xq'°) = (Xq ,(7 \) reads: 



D(k, E) 



G^{E) - P(k, E% - P(k, E)%f (X^) + t(k, E)\ " P(k, £0.(6) 



Here, Gq 1 ^) is the exact local Green function, Gq(E) = 5 a p/[E — (e n — £„')], 
S(A;, E) and P(k, E) are the self-energy and the strength operators, respectively. The 
presence of the strength operator is due to the redistribution of the spectral weight 
between the Hubbard subbands, that is an intrinsic feature of the strongly correlated 
electron systems. It is also should be stressed that S(/c, E) in Eq. ([6]) is the self-energy in 
the X-operators representation and therefore it is different from the self-energy entering 

— * 

the Dyson equation for the weak coupling perturbation theory for G\(k, E). 

To calculate the strength operator P(k, E) we use the zero-loop approximation given 
by the replacement: P a/3 (k,E) — > P a/3 = 5 a pF ai where F a ( n>n i) - 



X 



f 



is the occupation factor. Taking into account that in the considered paramagnetic 
phase (x^ ,a ^ = ^=^, ^X®' ^ = x, after all substitutions and treating all ^-independent 
terms as the chemical potential renormalization, the generalized Dyson equation for the 
Hamiltonian (TSl) becomes: 



D(k, E) 



(1 + x)/2 



E-(e - n) 



l+XJ. . 

2 k 



1+x t g 1-x 
2 E ct 2 



+ X{k,E) 



(7) 



To go beyond the Hubbard-I approximation we have to calculate S(fc, E). This was 
done in Ref. [27] using an equations of motion method for the X-operators [28]. The 
calculations resulted in the following expression: 



xN 



E 



X 



Jr. 



x 



E, 



1 + x 2tA 



k"1 



ct 



E, 



ct 



Kn 



u 



x 




X ^kH-q 



Here N is the number of vectors in momentum space, 
correlation function 




f,g 



2E' 



-irq 



and the kinematic correlation function 



K ir =J2e~ i{f ~ 9)ir (Xf X i 

f,9 



(8) 

the static spin-spin 
(9) 

(10) 
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Table 3. Doping-dependence of the effective t — t' — t" — J* model parameters (all 
values are in eV). Note, in the notations of Ref. [8] the hopping matrix elements tf g 
corresponds to —t°P, and i/ 9 = — t°f . Also, Jf g = 2tl /E ct . 



Parameter 



x =0.00 



0.05 



0.10 



0.15 



0.20 



0.30 




0.552 
0.054 
0.086 
0.463 
0.007 
0.012 
0.680 
0.085 
0.111 



0.560 
0.053 
0.087 
0.477 
0.007 
0.013 
0.691 
0.085 
0.112 



0.561 
0.050 
0.087 
0.484 
0.007 
0.013 
0.695 
0.081 
0.112 



0.572 
0.056 
0.070 
0.488 
0.007 
0.013 
0.699 
0.086 
0.113 



0.572 
0.054 
0.089 
0.492 
0.007 
0.013 
0.701 
0.084 
0.113 



0.567 
0.052 
0.088 
0.486 
0.007 
0.013 
0.697 
0.082 
0.112 



were introduced. 

Kinematic correlation functions ffTUj) are calculated straightforwardly using the 
Green function ([7j). The spin-spin correlation functions for the t — J model with 
the three-site correlated hoppings H 3 were calculated in Ref. [29J and the following 
expression for the Fourier transform of the spin-spin Green function was derived: 



where A^(uj) and magnetic excitations spectrum are given in Ref. [29] by the Eqs. (25) 

and (26), respectively 

The following results were obtained by the self-consistent calculation of the chemical 
potential fi, the spin-spin correlation functions ([9]) using Green function ( |TT|) . and the 
kinematic correlation functions (fTUj) using Green function ([7]) with the self-energy ([8]). 

Parameters of the effective t — t' — t" — J* model were obtained directly from 
the ab-initio parameters of the multiband model, Table [2J Their dependence on Ce 
concentration is presented in Table [31 Note, here we took Cu-4s orbitals into account 
implicitly through LDA+GTB method, as was described in Section However, within 
the Hubbard operators technique it is possible to take these orbitals into account 
explicitly using more cumbersome methods like the one introduced in Ref. [30J. 

Our results for the doping dependence of the kinematic and spin-spin correlation 
functions are shown in Fig. [31 Here, variable x' denotes the concentrations for which 
the one-electron parameters were calculated. Thus, x' = corresponds to the absence 
of the one-electron mechanism of doping dependence, while x' = x corresponds to 
the presence of this mechanism. However, the many-body mechanism is present in 
both cases. Note, the behavior of all correlation functions is almost identical for 
the cases of presence and absence of one-electron mechanism of doping dependence. 
Also note, the kinematic correlation functions, K n , possess a very nontrivial doping 
dependence. For low concentrations, x < 0.2, due to the strong magnetic correlations 
the hoppings to the nearest neighbors are suppressed leading to the small value of K\, 




(11) 
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Figure 3. (color online) Doping dependent evolution of the kinematic (upper panel) 
and the spin-spin (lower panel) correlation functions within the t — t' — t" — J* model. 
Index n enumerates the real space vectors connecting the neighboring sites: n = 1 for 
the nearest-neighbors, n = 2 for the next nearest neighbors, and so on. 



while K 2 and K 3 are not suppressed. Upon increase of the doping concentration above 
x ~ 0.2, magnetic correlations decrease considerably and the nearest-neighbor kinematic 
correlation function K\ increase, reviving the almost Fermi liquid behavior: K\ becomes 
largest of all K^s, while the magnetic correlation functions, C n , are strongly suppressed. 

So, we can clearly define one point of the crossover, namely, x m ~ 0.2. The system 
behavior is quite different on the different sides of these point, although there is no 
phase transition with symmetry breaking occurs. Apparently, this crossover is closely 
connected to the change of the Fermi surface topology with doping. Fermi surface 
evolution together with the quasiparticle dispersion is shown in Fig. HI For small x, the 
electron pockets around (±7r, 0) and (0, ±7r) points are present at the Fermi surface. 
Upon increase of the doping concentration these pockets become larger and merge 
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Figure 4. Doping-dcpcndcnt quasiparticle dispersion (on the left) and Fermi surface 
(on the right) in the spin-liquid phase of n-type cupratcs. The position of the chemical 
potential is denoted by the horizontal (green) line. The solid blue curves correspond 
to the calculations without taking the onc-clcctron mechanism of concentration 
dependence into account (x 1 — case) . The dashed red curves represent results of the 
calculations in which both many-body and one-electron mechanisms were considered 
(x' = x case). Note, the Fermi surfaces for both x' = and x = x cases are almost 
indistinguishable. 
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together at x = 0.2. For higher concentrations, the Fermi surface appears to be a 
large hole-like one, shrinking toward the (tt, tt) point. Thus, the topology of the Fermi 
surface changes at the same doping x m , where the point of crossover is situated. For 
the first time the "electronic transition" accompanying the change of the Fermi surface 
topology, or the so-called Lifshitz transition, was described in Ref. [31]. Now such 
transitions referred as a quantum phase transitions with a co-dimension= 1 (see e.g. 
paper [32]). 

Note, when the Fermi surface topology changes at a quantum critical concentration 
x m = 0.2 the density of states at the Fermi level also exhibit significant transformations. 
This results in the different behavior of the kinematic and magnetic correlation functions 
on the different sides of this crossover point. And, of course, the changes in the density 
of states at the Fermi level will also result in the significant changes of such observable 
physical quantities as the resistivity and the specific heat. 

Above the critical concentration, the Van Hove singularity in DOS is due to the 
flat dispersion around the (tt, 0) point. However, for x < x m this singularity is due to 
the states near of the (tt/2, tt/2) point. This is the result of the many-body interactions 
that can not be obtained within the one-electron band theory. 

Concerning the role of the short range magnetic order and three-site hopping terms 
in the n-type cuprates we would like to stress that due to the scattering on the magnetic 
excitations with the AFM wave vector Q = (tt, tt) the states near the (tt, it) point are 
pushed above the Fermi level, and the local symmetry around the (tt/2, tt/2) points is 
restored for low doping concentrations (see Fig. H|). In other words, the short range 
magnetic order "tries" to restore the AFM symmetry around (tt/2, tt/2) point. In our 
calculations, the short range magnetic fluctuations are taken into account via the spin- 
spin correlation functions (Q. 

Now we proceed to the comparison of the one-electron and many-body mechanisms 
of the doping dependence. In Fig. @] the quasiparticle dispersion without one-electron 
mechanism is shown by the solid blue curves. Apparently, its difference from the case 
when both one-electron and many-body mechanisms are present is negligibly small. In 
the latter case the bandwidth is slightly renormalized while band dispersion retains the 
same character. Moreover, the Fermi surfaces for both cases are very similar and the 
quantum phase transition will be at the same concentration, x m = 0.2. 

6. Conclusion 

In the present work we report the combined investigation of the one-electron and the 
many-body mechanisms of the electronic structure doping dependence for the High-T c 
compound Nd2- x Ce 2; Cu04. The electronic structure calculations were performed within 
the hybrid LDA+GTB scheme. For the both antiferromagnetic and paramagnetic spin- 
liquid phases we demonstrate that the main effect on the electronic structure is provided 
by the many-body mechanism, whereas the one-electron contribution leads to a rather 
small quantitative modifications which do not change the picture qualitatively. The 
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role of the many-body mechanism is very important because of the strong electronic 
correlations present in the underdoped cuprates. 

Acknowledgments 

This work was supported in part by RFBR Grants 06-02-16100 (M.M.K., V.A.G., 
S.G.O.), 06-02-90537-BNTS (M.M.K., V.A.G., S.G.O., I.A.N., Z.V.P.), 05-02-16301 
(I.N.), 05-02-17244 (I.N.), by the joint UrO-SO project 74, and programs of the 
Presidium of the Russian Academy of Sciences (RAS) "Quantum macrophysics" 
and of the Division of Physical Sciences of RAS "Strongly correlated electrons 
in semiconductors, metals, superconductors and magnetic materials". M.M.K. 
acknowledge support form INTAS (YS Grant 05-109-4891). I.N. and Z.P. acknowledge 
support from the Dynasty Foundation and International Center for Fundamental Physics 
in Moscow program for young scientists. I.N. appreciate the support from the grant of 
President of Russian Federation for young PhD MK-2242. 2007.02. Z.V. is supported by 
UrO grant for young scientists. 

References 

[1] Massidda S, Hamada N, Yu Jaejun, Freeman A J 1989 Physica C 157 571 

[2] Yu Jaejun, Freeman A J 1991 Physica C 173 274 

[3] Yu Jaejun, Freeman A J 1991 J. Phys. Chem. Solids 52 1351 

[4] Mishonov T, Pcncv E 2000 J. Phys.: Condens. Matter 12 143 

[5] Andersen O K, Liechtenstein A I, Jepsen O, Paulsen F 1995 J. Phys. Chem. Solids 56 1573 
[6] Andersen O K, Savrasov S Y, Jepsen O, Liechtenstein A I 1996 J. Low Temp. Phys. 105 285 
[7] Pavarini E, Dasgupta I, Saha-Dasgupta T, Jepsen O, Andersen O K 2001 Phys. Rev. Lett. 87 
047003 

[8] Korshunov M M, Gavrichkov V A, Ovchinnikov S G, Nekrasov I A, Pchclkina Z V, Anisimov V I 

2005 Phys. Rev. B 72 165104 
[9] Gavrichkov V A, Ovchinnikov S G, Borisov A A, Goryachcv E G 2000 Zh. Eksp. Teor. Fiz. 118 

422 [JETP 91 369] 

[10] Gavrichkov V A, Ovchinnikov S G 2004 Zh. Eksp. Teor. Fiz. 125 630 [JETP 98 556] 

[11] Kamiyama T, Izumi F. Takahashi H, Jorgensen J D, Dabrowski B, Hitterman R L, Hinks D G, 

Shaked H, Mason T O, Seabaugh M 1994 Physica C 229 377 
[12] Paulus E F, Ychia I, Fucss H, Rodriguez J, Vogt T, Strobel J, Klauda M, Saemann-Ischenko G 

1990 Solid State Commun. 73 791 
[13] Andersen O K, Jepsen O 1984 Phys. Rev. Lett. 53 2571 
[14] Andersen O K, Pawlowska Z, Jepsen O 1986 Phys. Rev. B 34 5253 
[15] Singh D J 1991 Phys. Rev. B 44 7451 
[16] Emery V J 1987 Phys. Rev. Lett. 58 2794 

[17] Varma C M, Smitt-Rink D. Abrahams E 1987 Solid State Commun. 62 681 
[18] Andersen O K, Saha-Dasgupta T 2000 Phys. Rev. B 62 16219(R) 

[19] Anisimov V I, Korotin M A, Nekrasov I A, Pchelkina Z V, Sorclla S 2002 Phys. Rev. B 66 100502 
[20] Gunnarsson O, Andersen O K, Jepsen O, Zaancn J 1989 Phys. Rev. B 39 1708 
[21] Hubbard J C 1963 Proc. R. Soc. A 276 238 
[22] Zaitsev R O 1975 Sov. Phys. JETP 41 100 

[23] Izumov Yu, Lctfullov B M 1991 J. Phys.: Condens. Matter 3 5373 



Dominance of many-body effects over one-electron mechanism in Nd2- x Ce x Cu04 15 

[24] Ovchinnikov S G, Val'kov V V 2004 Hubbard Operators in the Theory of Strongly Correlated 

Electrons (Imperial College Press, London-Singapore) 
[25] Hubbard J C 1964 Proc. R. Soc. A 277 237 

[26] Ovchinnikov S G, Borisov A A. Gavrichkov V A, Korshunov M M 2004 J. Phys.: Condens. Matter 
16 L93 

[27] Korshunov M M, Ovchinnikov S G 2006 Eur. Phys. J. B 57 271 

[28] Plakida N M, Yushankhai V Yu, Stasyuk I V 1989 Physica C 162-164 787 

[29] Val'kov V V, Dzcbisashvili D M 2005 Zh. Eksp. Teor. Fiz. 127 686 [JETP 100 608] 

[30] Digor D F, Entcl P, Moskalenko V A, Plakida N M 2006 Theor. and Math. Phys. 149 1382 

[31] Lifshitz I M 1960 Sov. Phys. JETP 11 1130 

[32] Volovik G E 2006 Acta Phys. Slov. 56 49; cond-mat/ 0601372 

[33] Plakida N M, Oudovcnko V S 2007 Physica C 460-462 993 

[34] Plakida N M, Oudovcnko V S 2007 Zh. Eksp. Teor. Fiz. 131 259 [JETP 104 230] 



